Two-dimensional scattering and bound states of polar molecules in bilayers 
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Low-energy two-dimensional scattering is particularly sensitive to the existence and the properties 
of weakly bound states. We show that interaction potentials V(r) with vanishing zero-momentum 
Born approximation J drrV(r) — lead to an anomalously weak bound state which crucially mod- 
ifies the two-dimensional scattering properties. This anomalous case is especially relevant in the 
context of polar molecules in bilayer arrangements. 
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I. INTRODUCTION 

Ultracold atomic gases are many-body systems, but 
most of their fundamental properties originate from the 
underlying two-body problem, given by an interaction 
potential which is typically considered as short-range. 
The situation is completely different in dipolar gases due 
to the long-range character of the dipole-dipole inter- 
action 0, [2j]. Polar molecules with a potentially large 
electric dipole moment constitute a particularly excit- 
ing dipolar gas. However, achieving quantum degener- 
acy is handicapped by exothermic chemical reactions Q . 
The latter may be avoided by confining the gas in two- 
dimensional geometries if the dipoles are polarized per- 
pendicular to the trap plane, due to the repulsive char- 
acter of the dipole-dipole interaction [J]. 

Bilayer (and in general multi-layer) arrangements of 
polar molecules offer the possibility of stability against 
inelastic reactions and give rise to inter-layer pairing due 
the dipole-dipole force [5|-|7[ . A dipole in layer 1 interacts 
with a dipole in layer 2, both dipoles oriented perpendic- 
ularly to the layers, by the potential 



V dd (r) = U (r 2 - 2)/(r 2 + if' 2 . 



(1) 



Here r is the relative in-plane distance between the two 
dipoles, the inter- layer distance A is set to 1, and Uq 
is a positive dimensionless coupling constant. We have 
Ua = md 2 /h 2 X and the unit of energy is Eq = h 2 /mX 2 , 
with m the mass and d dipole moment of the molecule. 
This potential is attractive at short distances and repul- 
sive at large distances, fulfilling the peculiar condition 
/ drrVdd(r) — 0, i.e. its zero-momentum Born approx- 
imation vanishes 0-@- Inter- layer interactions of a dif- 
ferent type may be attained in binary mixtures, where 
one of the species is confined in a bilayer while the other 
moves freely. The free species mediates a 2D interaction 
with a RKKY-like potential Q3, which may as well have 
a vanishing zero-momentum Born approximation. 

Low energy 2D scattering, which determines the prop- 
erties of 2D quantum gases pH - [l3j . is particularly sen- 
sitive to the existence and properties of weakly bound 
states. Although 2D scattering |l#tl§ ] and weakly bound 
states [l9l - |2lj have been intensively studied, little is 



known for the case when J drrV{r) — 0. The binding 
energy for weakly coupled bound states in this case was 
calculated in Ref [l|. However, a detailed investigation 
of the binding energy at larger coupling and of the low- 
energy scattering properties is still lacking. 

In this Brief Report we discuss the low-energy scatter- 
ing and weakly bound states for radial potentials V(r) 
satisfying J drrV(r) = 0. We extend (using an alter- 
native method) the expression derived in Ref. for 
the binding energy of the weakly bound state. We show 
that the presence of this anomalously weak bound state 
modifies significantly the scattering amplitude compared 
to the usual case of potentials with non- vanishing zero- 
momentum Born approximation. As an example, we spe- 
cialize for the potential Vdd appearing in bilayer gases of 
polar molecules and check the validity of the obtained an- 
alytical expressions using exact numerical calculations. 

This paper is organized as follows. In Sec. UH we in- 
troduce the Jost function formalism to study the binding 
energy of weakly bound states and evaluate the general 
expressions for Vdd- In Sec. IIIII we discuss the modifica- 
tions introduced in the two-dimensional scattering prop- 
erties if J drrV{r) — and study this for Vdd in detail. 
Our conclusions are summarized in Sec. IIVI 



II. WEAKLY BOUND STATES 

Ref. [19( studied the bound states of the Schrodinger 
equation {—V 2 + V(r)}ip — etp for potentials of the 
form V(r) = Uov(r), where as above Uq denotes a pos- 
itive dimensionless coupling constant characterizing the 
potential strength [26]. For the case of weak coupling 
(Uq —> 0) it was shown that a shallow bound state al- 
ways exists if J d 2 rv(r) < 0, but there is no bound state 
if J d 2 rv(r) > 0. Furthermore, it was shown that for 
/ d 2 rv(r) < the binding energy of the shallow bound 
state is of the form eb ~ — exp[47r/ J d 2 rV(r)], as ex- 
pected from e.g. Ref. 22[. However, for our case of in- 
terest / d 2 rv(r) — 0, the binding energy acquires the 
anomalous form ~ — exp[l/cC/g], with 



1 

8^- 



c = — I d 2 r 



J d 2 r'v{r) ln|r-r>(r'). 



(2) 



2 



It can be shown that c < for any v(r). For the potential 
Vdd(^) appearing in a bilayer system of polar molecules, 
one obtains c = —1/8, and hence the binding energy 
becomes: 



Ad 



exp{-8/U 2 } 



(3) 



However, a numerical calculation of the binding energy 
(see below) shows that this result is not very accurate 
even for very small Uq. This motivates us to derive a 
more accurate analytic expression for the binding energy 
which remains valid for larger Uq- 

The two-body scattering problem for a radially sym- 
metric potential V(r) in two dimensions is described by 
the Schrodinger equation 



ar z r dr 



where all quantities are dimensionless and 4>{r) is the ra- 
dial wavefunction. Only s-waves are considered since we 
are only interested in low-energy properties. Following 
Ref . [23| , we employ the Jost function formalism to study 
the scattering problem and the shallow bound states. 
The definition of the Jost function T(k) in the 2D case 
may be found in Ref. [23|. The properties of T(k) are 
similar to those of the Jost function for the 3D case [24| . 
The scattering phase shift 5(k) is related to J-(k) by 



tan<5(fc) = — 



Im T{k) 
Re F{k) ' 



the scattering amplitude is 



/(*) 



tan<5(fc) 
1 — ttan(5(fc) 



(5) 



(6) 



and a = (4/k)\f(k)\ 2 = (4/fc) sin 2 5(k) is the total 2D s- 
wave cross-section. For complex k, the zeros of the Jost 
function on the positive imaginary axis, ^{ia) = with 
a > 0, are the bound states of the potential with binding 
energy e b = -a 2 . 

The following integral representation of the Jost func- 
tion -F(fc) will be employed below: 



T(k) = l + e^/^—J^ dr^V(r)J (kr)f (kr), (7) 
where fa(k,r) satisfies the integral equation 



iirkr 



H (kr) + dsg(k, r, s)V(s)f (k, s), (8) 



fo(k,r) 
with 



g(k,r,r') = -V^ [J {kr)Y (kr') - J {kr')Y (kr)\ . 

In the previous expressions Jo, Y are Bessel functions, 
and Hq is the Hankel function of the first kind. 
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FIG. 1: (Color online) Binding energy for Vdd as a function 
of the dipole strength t/o calculated numerically (dots), from 
equations (|16[) (solid) and ((3| (dashed). 

In the following we are interested in determining the 
weakly bound states. To this aim we expand the Jost 
function for small k: 



F{k) = Alnk + B -i-A, 



(9) 



where A, B are real constants independent of k. These 
constants can be represented by infinite series of the form 

A = E A ii B = 1 ~ E b j E3 ■ The nrst terms of the 
Aj series are of the form: 



drrV{r) 



(10) 



drrV(r) dssV(s) In (^-J , (11) 



A x 
A 2 



whereas those of the Bj series are: 

/•oo 

Bi= drrV(r) ln(r) - C\ (12) 

B 2 = J drrV(r) J dssV(s) In ^- J ln(s) - C 2 , (13) 

with Ci — Ai ln(e 7 /2) and 7 ps 0.577 the Euler constant. 

As mentioned above, the bound states are given by 
the zeros of T{k) on the positive imaginary axis. Using 
(j9]) with k — ia we hence obtain the expression of the 
binding energy: 



exp 



4 



(14) 



Note that Eq. (ITil) is valid as long as the binding energy 
is small enough such that the logarithmic term dominates 
the Jost function. For potentials with J drrV(r) < 0, it 
is sufficient to take A m A\ and B « 1, recovering the 
expression e& ~ exp(4-7r/ J d 2 rV(r)) for small Uq. 

However, for the case of potentials such that 
J drrV(r) — 0, we have A\ — 0, and the first non- 
vanishing term is A A 2 , B 1, providing an alter- 
native derivation of Eq. A more precise formula is 
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FIG. 2: (Color online) Scattering phase shift for Vdd as a 
function of k for different Uo calculated numerically (dots) 
compared to the log-behavior (|20[1 for Uo = 1.0 (dash-dotted, 
green), to formula (|17[1 for Uo = 0.6 (solid, orange) and to the 
second Born approximation (|21[) for Uo = 0.6 and Uo — 0.2 
(dashed, cyan). 



obtained by including higher-order terms: 



exp < —2 



l-Bx-B 2 - 
A 1 + A 2 + . . 



(15) 



For the case of the inter-layer dipole-dipole potential 
Vdd) the integrals can be carried out analytically to find 
a corrected expression for the binding energy 



rid. 



expi-^ 



1-Uo 



(16) 



Figure Q] compares the numerical result for the bind- 
ing energy for Vdd (obtained directly from the 2D 
Schrodinger equation) with the analytical expressions of 
Eq. Q and ([TB"]) . Note that whereas Eq. ([3]) provides a 
relative inaccurate approximation even at rather low Uq, 
the newly derived expression (1161) is in excellent agree- 
ment with the numerics, all the way to Uq < 1.2. 

Finally, we note that for large Uq the binding energy 
for Vdd can be determined by a variational calculation, 



giving 25] e 



dd 



-2U + 4^3LV2 - 15/4 , which 



cides with the numerics only for Uq > 5. 



III. SCATTERING PHASE SHIFT 

We have shown above that the Jost function formal- 
ism is particularly useful for the analysis of weakly bound 
states. In this section we employ this formalism for the 
study of 2D scattering and in particular for the calcula- 
tion of the s-wave scattering phase shift 8(k). 

An approximate expression of the Jost function T{k) 
for small Uq is obtained by iterating twice the integral 
equation ([5]). Note that we keep all orders in k. The 
resulting scattering phase shift follows from the rela- 
tion ©: 



tan<5(fc) ; 



f!j,/(fc)-T {hj,JY{k)-lj Y ,jj{k)) 



-I 



JY 



(I 



JJ,YY 



(*)- 



(17) 




FIG. 3: (Color online) Qualitative sketch of the regimes of k 
and Uo where the scattering phase shift for Vdd can be approx- 
imated by the logarithm Eq. (|20[) (green), the Born approxi- 
mation Eq. Q22p (yellow) , and the second Born approximation 
Eq. (2D (orange). 



where we have introduced the notation 

f-OC 

Ifg= drrV(r)F(r)G(r), 
Jo 

I 



(18) 



i i; P q=I drrV(r)F(r)G(r)l dssV(s)P(s)Q(s), (19) 



and J, Y stand for Jo(kr) and Y (kr). 

For small k it is possible to simplify Eq. (TTT)) . Employ- 
ing the logarithmic expression (9]), the relation (SJ, and 
the expression for the binding energy (|14[) , we recover the 
well known logarithmic expression (see e.g. [l2l ]) 



tan 5(k) 



1 , k 2 
-In — 

7t \e b \ 



(20) 



characteristic of 2D scattering, which relates the scat- 
tering shift and the binding energy of the weakly bound 
state. However, for the case J drrV{r) = the binding 
energy |et,| can become anomalously small, and hence the 
expression at the right-hand side of Eq. (|20|) can become 
very small for reasonable k. In this case, it is not any 
more the leading term for the low-energy scattering. 

On the other hand, for large enough k and small Uq, 
the first integral Ijj(fc) in Eq. (|T7)) dominates, and we 
recover, as expected for sufficiently large k, the Born 
approximation. Therefore, formula (1171) interpolates 
smoothly between the correct low-energy and the cor- 
rect high-energy behavior. This suggests that it may be 
valid, at least qualitatively, even for large Uq. By ex- 
panding Eq. (|17p in powers of Uq, one recovers the Born 
series, which reads to second order: 



tanJ(fc) 



j j 



IJJ,JY- 



(21) 



It is interesting to discuss these results for the case 
of the interlayer dipole-dipole potential Vdd- We have 
computed the scattering phase-shift numerically from the 
Schrodinger equation and compared it to the results of 
Eq. (fi~7) . obtaining that Eq. (fT7| provides the correct 



4 



scattering phase-shift with excellent accuracy, at least in 
the range 0.03 < k < 5 and 0.05 < Uq < 2.0. Hence 
Eq. ([P7)l is a good approximation not only for t/ « 1, 
but also for Uq ~ 1. Of course, if the interactions are too 
large (Uq 3> 1), the contributions from higher iterations 
of the integral equation ((8]) become more important and 
Eq. (fT7|) looses its accuracy. 

Figure[2]compares the numerical results for the scatter- 
ing phase-shift for Vdd with the limiting cases provided by 
expressions (|20[) and (|2"Tj) for small k and different values 
of Uq. It is seen that for Uq ~ 1, the scattering phase-shift 
is best approximated by the logarithmic expression (|20p , 
and for Uq <C 1 by the second Born approximation (f2Tj) . 
For intermediate values of Uq none of the limiting cases 
is accurate and the full expression (|17[) must be used. 
We sketch in Fig. [3] qualitatively the regimes of k and 
Uq, where the logarithm (|20p. the first Born approxima- 
tion (|22[) . and the second Born approximation (f2"TT) arc 
good approximations, as obtained by comparison with 
the numerical solution. Note that, excluding unreason- 
ably small k, the logarithmic form (|2T))) is just valid for 
k < 1 and the window 0.7 < Uq < 2.0. 

Finally, we note that the first Born approximation for 
Vdd can be evaluated exactly analytically: 



tan,5 dd (fc) 



;U 



Ik 



•2fc(Li(2fc) -h{2k)) 



(22) 



IV. CONCLUSIONS 

In conclusion, two-dimensional radial interaction po- 
tentials V(r) with a vanishing zero-momentum Born ap- 
proximation, J drrV(r) = 0, result in interesting physics 
crucially different from purely attractive or purely repul- 
sive potentials. Using the Jost function formalism, we 
have derived an expression for the binding energy as a 
function of the potential strength Uq, which remains ac- 
curate for a wide regime of Uq values. Moreover, we have 
investigated the scattering amplitude in different param- 
eter regimes. In particular, we have shown a significant 
deviation of the scattering behavior in comparison with 
potentials with JdrrV(r) ^ 0, due to the anomalously 
low binding energy of the weakly bound state. 

These results are of particular importance in the 
physics of two-dimensional systems, and more specifically 
on two-dimensional ultracold gases. Standard theories, 
in particular the theory of BCS-BEC crossover [l2l. [l3|. 
are based on the fact that the scattering amplitude pos- 
sesses the logarithmic dependence ([2"U)) . These results 
are therefore modified if the potential has a vanishing 
zero-momentum Born approximation. This has particu- 
larly important consequences for the properties of a gas 
of polar Fermi molecules confined in a bilayer geometry, 
including inter-layer pairing [7J. 



where Li is the modified Struve function. The second 
Born approximation (|21j) can be expanded for small k 

tan(5 dd O) 2U k - TrU a k 2 + ±[/ 2 ± . . . , (23) 

which gives a maximum in the scattering phase shift at 
k sa 1/tt as observed in the numerical results. 
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